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Abstract 

The past three decades of investigation on nuclear physics and pulsar astro- 
physics have seen gradual recognition that elasto dynamic approach to the contin- 
uum mechanics of nuclear matter provides proper account of macroscopic motions 
of degenerate Fermi-matter constituting interior of the nuclear material objects, the 
densest of all known today. This paper focuses on one theoretical issue of this de- 
velopment which is concerned with oscillatory behavior of a viscoelastic solid globe 
in the regime of quasistatic, force-free, non-compressional oscillations less inves- 
tigated in the literature compared to oscillations in the regime of standing shear 
waves. We show that in this case the problem of computing frequency and lifetime 
of spheroidal and torsional modes of non-radial shear vibrations damped by viscos- 
ity can be unambiguously resolved by working from the energy balance equation 
and taking advantage of the Rayleigh's variational method. The efficiency of this 
method is demonstrated by solid globe models of nuclear physics and pulsar astro- 
physics dealing with oscillations of a spherical mass of a viscoelastic Fermi-solid with 
homogeneous and non-homogeneous profiles of the bulk density, the shear modulus, 
and the shear viscosity. 
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1 Introduction 



An understanding dynamical laws governing macroscopic motions of degenerate nucle- 
onic material - the continuum mechanics of nuclear matter - is important for developing 
interconnected view of the nuclear physics and pulsar astrophysics. The nearly linear 
dependence of the nucleus volume and the nucleus binding energy upon number of nu- 
cleons, meaning the saturation of internucleon forces, unambiguously indicate that this 
"most conspicuous feature of nuclei" (Bethe 1999) is the same as it is for condensed matter. 
The condensed matter, as is known, comes in two forms - liquid and solid. Accordingly, 
the question which of two fundamental models of material continua - the fluid-mechanical 
(relying on equations of hydrodynamics highlighting fluidity as the major dynamical prop- 
erty of liquid state of condensed matter) or the solid-mechanical (resting on equations of 
elastodynamics accentuating reversal elastic deformations as basic dynamical property 
of solid state of condensed matter) can provide adequate description of the observable 
behavior of the nuclear matter objects has been and still is central to the program on the 
study of material properties of superdense substance constituting interior of atomic nuclei 
and neutron stars. 

It is of common knowledge today that the earlier liquid drop model of nucleus, implying 
the fluid-mechanical pattern of continuous nuclear matter and complete disorder in the 
behavior of constituting particles, has fallen in disfavour after discovery of structural order 
in proton-neutron composition of nucleus which received theoretical interpretation in the 
nuclear shell model. This microscopic signature of ordered clusterization of nucleons 
together with such macroscopic one as the deformed equilibrium shapes indicate that 
the atomic nucleus bears more resemblance to an ultrafine particle of a viscoelastic solid 
than a drop of liquid which is characterized by spherical equilibrium shape. Such an 
understanding is clearly traced in the macroscopic approaches to the nuclear collective 
dynamics and asteroseismology of neutron stars developed in the past three decades that 
have been dominated by mathematical methods and physical notions of the classical 
elastodynamics testifying gradual receding of the fluid-mechanical paradigm of the nuclear 
continuum mechanics into the solid mechanical one. While the microscopic nature of 
structural order underlying substantially quantum clusterization of nucleons in the nuclear 
matter objects, which has definitely nothing to do with classical picture of lattice-like 
crystallization of atoms in ordinary solids, is still to be clarifiecQ, an independent inference 
of current nuclear physics and pulsar astrophysics is that the elastodynamical model of 

In current investigations on quantum condensed media, the term Fermi-solid (which is, as far as we 
know, due to Guyer (1973)) is applied to accentuate the difference between single-particle behavior of 
microconstituents in quantum and ordinary solids (e.g. Polturak and Gov 2003). 
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nuclear continuum mechanics, lying at the base of the solid globe model for an atomic 
nucleus and the solid star model for a pulsar, provides proper macroscopic account of 
data on giant-resonant oscillatory response of nuclei and non-radial pulsations of neutron 
stars. 

By now there are quite cogent theoretical arguments showing that the asteroseismology 
of neutron stars, studying restless behavior of pulsars and magnetars, can be properly 
understood working from the model of a solid staj^l, rather than the liquid star model 
underlying current understanding of seismic activity of the main sequence stars. Also, it 
is believed the much more denser strange matter of quark stars (e.g. Weber 2005) is in 
the solid state too (Xu 2005; Owen 2005). In broad terms, by a solid star is understood 
a self-gravitating mass of viscoelastic solid continuum, highly robust to compressional 
deformations, whose motions are governed by the laws of Newtonian elastodynamics. In 
pulsar astrophysics, a model of a solid star undergoing global and surface (crustal) elastic 
vibrations has been invoked to explain microspikes detected in the windows of the rotation 
driven main pulse train of the pulsar radio emission (e.g. Van Horn 1980, McDermott, 
Van Horn and Hansen 1988; Strohmayer, Cordes and van Horn 1993; Bastrukov, Weber 
and Podgainy 1999). The core-crust model of a solid star was found to provide proper 
account of the neutron star quakes manifested by glitches in the regular pulsating emission 
of radio pulsars (e.g. Pines, Shaham and Ruderman 1974; Canuto and Chitre 1974; 
Franco, Link and Epstein 2000) and gamma bursts of pulsars and magnetars (Blaes, 
Blandford, Goldreich, Madau 1989; Cheng, Epstein, Guyer and Young 1996). In this 
model the neutron star interior is thought of as a stiff outer crust covering a denser core 
which are separated by thin inner crust whose matter operates like a lubricant facilitating 
tangential shear displacements of the outer crust against core. This picture suggests 
that the post-quake recovering of glitching radio pulsars and gamma-bursting magnetars 
to the regime of quiescent pulsed radiation is accompanied by the torsional, that is, 
differentially-rotational elastic vibrations of the neutron star crust against core. Over 
the years the torsional pulsations of neutron stars have gained considerable attention in 
different aspects of radiative activity of pulsars and magnetars (e.g. Ruderman 1969; Van 
Horn 1980, Hansen and Chioffi 1980; Bastrukov, Podgainy, Yang, Weber 2002). This 
kind of vibrations, generic to oscillatory behavior of an elastic solid globe, is one of the 
topical issues in current discussion of periodicity in bursting activity of X-ray pulsars 

2 To the best of our knowledge the term "Solid Stars" for neutron stars has been suggested by M. 
Ruderman in the infancy of pulsar astrophysics, as pointed out in (Harwit 2006). A current view on the 
neutron star structure is extensively discussed in (Shapiro and Teukolsky 1983; Weber 1999; Glendenning 
2000; Lattimer and Prakash 2001). 
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(e.g. Lin and Chang 2005, Strohmayer and Watts 2006; Strohmayer and Bildsten 2006; 
Bisnovatyi-Kogan 2007). 

In nuclear physics, the signatures typical of oscillatory behavior of viscoelastic solid 
globe have been discovered in the high-energy collective response of atomic nucleus by 
giant resonances (Bertsch 1974) - the fundamental modes of vibrational excitations generic 
to all stable nuclei of the periodic table (e.g. Harakeh and van der Woude 2003; Richter 
2005). Theoretical treatment of these modes in terms of spheroidal (electric resonances) 
and torsional (magnetic resonances) elastic vibrations of a solid sphere has been and still 
is the subject of extensive research by different methods of the nuclear collective dynamics 
(e.g. Semenko 1977; Sagawa and Holzwarth 1978; Nix and Sierk 1980; Wong and Azziz 
1981; Hasse, Ghosh, Winter and Lumbroso 1982; Stringari 1983; Balbutsev and Mikhailov 
1988; Kolomietz 1990; Di Toro 1991; Bastrukov, Misjcu and Sushkov 1993; S. Misjcu 2006 
and references therein). One of the main goals of these investigations is to get insight into 
material laws, that is the rheology of the nuclear matter. The torsional oscillatory mode 
is unique to the nuclear solid globe model relying on equations of the solid mechanics. 
Thereby, the macroscopic interpretation of the giant magnetic resonances discovered in 
the photoelectron experiments (Raman, Fagg and Hicks, 1991; Luttge, von Neumann- 
Cosel, Neumeyer and Richter 1996) as manifestation of the nuclear twist (Holzwarth 
and Eckart 1977; Bastrukov and Gudkov 1992), that is, torsional differentially-rotational 
elastic oscillations of a nucleus (Bastrukov, Molodtsova and Shilov 1995; Bastrukov, Libert 
and Molodtsova 1997) serves as the most strong argument in favor of the elastodynamical 
model of nuclear continuum mechanics. The equations of fluid mechanics underlying 
canonical liquid drop model of nucleus (e.g. Greiner and Maruhn 1996) have no solutions 
describing nuclear response by differentially rotational oscillations and this is the reason 
why the liquid drop model excludes vibrational states of magnetic type, as was pointed 
out long ago in (Rowe 1970). 

This paper focuses on one subtle issue of general solid-mechanical treatment of os- 
cillatory behavior of the nuclear matter objects in terms of vibrations of a viscoelastic 
solid globe that has been called into question by the above development of both the 
nuclear physics and pulsar astrophysics. Specifically, it is concerned with the regime of 
quasistatic, force-free, oscillations less investigated in the literature as compared to the 
solid globe oscillations in the regime of standing shear waves. The most striking feature 
of the quasistatic regime is that the restoring force of Hookean elastic (reversal) stress 
and dissipative force of Newtonian viscous (irreversal) stress entering the basic equation 
of solid mechanics turn to zero (from what the term quasistatic oscillations is derived), 
but the material stresses themselves and the work done by these stresses in the bulk of an 
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oscillating solid globe do not. Based on this observation we show that in the case of force- 
free fluctuations of stresses the frequency and lifetime of both spheroidal and torsional 
modes in a viscoelastic solid globe can be computed by taking advantage of the energy 
variational principle relying on the equation of energy conservation, not on the equation 
of elastodynamics as is the case of oscillations in the regime of standing waves. 

The paper is organized as follows. In Section 2 a brief outline is given of solid- 
mechanical equations modeling oscillatory behavior of a viscoelastic solid with focus on 
the problem of oscillatory response of a perfectly elastic (non-viscous) sphere in the regime 
of standing shear waves that dates back to the canonical work of Lamb (1882). Section 
3 is devoted to the solution of the problem of quasistatic oscillations of a viscoelastic 
solid globe by Rayleigh's variational method and its application to the nuclear giant 
resonances. The efficiency of this method is demonstrated in Section 4 considering several 
mathematical models of astrophysical interest. The practical usefulness of these models 
is briefly discussed in Summary and details of calculations are presented in Appendices. 

2 Elastodynamics of the solid globe oscillations in the 
compression free regime of standing shear waves 

In approaching the topic like this, it seems worthy to start with a brief outline of physi- 
cal ideas and mathematical equations underlying the standard continuum model of solid 
mechanics or elastodynamics (e.g. Love 1944; Morse and Feshbach 1953; Eringen and 
Suhubi 1975; Landau, Lifshits, Kosevich and Pitaevskii 1995). The basic dynamical vari- 
able characterizing the state of motion of solid continuous medium is the field of material 
displacement Ui(r,t) whose emergence in the volume of solid object is associated with its 
response to perturbation of mechanically equilibrium state in which Ui = 0. The canonical 
continuum model of an isotropic viscoelastic solid rests on two fundamental, empirically 
verified, statements. First is the generalized Hookean law of elastic deformations which 
is expressed by linear relation between the symmetric tensors of elastic (reversal) stress 
Pik and shear strain or deformation u^: 



ft. = + (« - |m) UjA* = 2, ( u, - ^» ) + 

Uik = -(VjMfc + VfcWj) U kk = VfcWfc. (2) 

Hereafter by the repeated indices is denoted summation over all values. The trace of the 
stain tensor, is the measure of fraction change in the volume of small elementary 
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mass of solid continuum u kk = 5V/V = —5p/p and hence 

bp = -pu kk u kk = V k u k . (3) 

The second fundamental statement of elastodynamics is the Newtonian law of viscosity 
which is described by similar tensorial form of constitutive equation expressing linear re- 
lation between the viscous stress tensor 7i ik and the rate-of-strain tensor u ik (e.g. Landau, 
Lifshits, Kosevich and Pitaevskii 1995) 

n ik = 2r]ii ik + ( C - \~n J u u 5 ik = 2rj \ u ik - \ujj8 ik \ + ( Ujj5 ik (4) 



3 7 ' V 3 

Uik = ^(VjMfc + VfeMj) U kk = VfcUfc. (5) 

In the analysis of oscillatory behavior of a viscoelastic solid, the strength's characteristics 
of viscoelastic deformations - the shear modulus p, the bulk modulus k (heaving, both p 
and k, physical dimension identical to that for the pressure P), the shear viscosity r\ and 
the bulk viscosity ( (in general local functions of position) are regarded, together with the 
bulk density p, as input parameters. These are specific to each given sort of viscoelastic 
material whose values are measured in experiments on propagation of material waves (so 
called primary or P-waves and secondary or S*-waves) in an isotropic solid. The second 
law of Newtonian dynamics for an isotropic viscoelastic medium is expressed by the Lame- 
Navier equation 

pUi = V k p ik + Vfc7T ifc . (6) 

Henceforth the overdot denotes partial differentiation with respect to time. The energy 
conservation in the process of viscoelastic distortions is controlled by the equation 

d pv? 

~Q t — = ~\Pik + Kik] U ik + Vfc[(p ifc + TTifc) Ui] (7) 

which results from after multiplication with Wj. 

From above it follows that non-compressional response of a viscoelastic solid to short- 
time applied load, which is not accompanied by local fluctuations in density, 5p = 
—pukk = 0, the tensors of viscoelastic stresses are reduced to 

p ik = 2pu ik n ik = 2r]u ik u kk = V k u k = 0. (8) 

When studying the fast mechanical response of a viscoelastic solid by waves of mate- 
rial displacements it is customarily assumed that intrinsic distortions are dominated by 
gradients in the tensor of strains and rate-of-strains 

\u ik V k p\ « \pV k u ik \ \u ik V k f]\ « \r]V k u ik \ (9) 
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implying 



ff = VkPik = 2u ik V k p + 2/iVfc u ik fiV 2 Ui 
f t N = V k Tc ik = 2u ik V k n + 2nW k u ik « nW 2 iii. 



(10) 
(11) 



so that parameters of viscoelasticity are regarded as constants. In this approximation the 
basic equation of solid mechanics is reduced to pui = pV 2 Ui + nV 2 Ui. 

2.1 Shear oscillations of a perfectly elastic solid globe in the 
regime of standing waves of material displacements 

To accentuate the subject of our investigation, we pause for a while on the case of perfectly 
elastic (n = 0) solid globe with uniform profile of shear modulus p, = constant. In this 
case equation of elastodynamics is reduced to the canonical wave equation pu = p V 2 u 
showing that perfectly elastic solid can transmit compressionless plane-wave perturbation 
u = u exp i(uot — k r) by transverse wave (V ■ u = — > k ■ u = 0) of elastic shear traveling 
with the speed cf = p/p and characterized by the dispersion free law of propagation 
uj 2 = c 2 k 2 . This suggests the problem of oscillatory response of a perfectly elastic solid 
globe can be solved by working from Helmholtz equation for the standing shear waves 
V 2 u + kfu = 0. As was mentioned, in this form the problem of oscillatory response 
of a perfectly elastic sphere has been posed and solved by Lamb (1882). Specifically, 
it has been suggested that oscillatory modes in an elastic sphere, termed as spheroidal 
(s — mode) and torsional or toroidal (t — mode), can be unambiguously classified by two 
general, regular in origin, solutions to the vector Helmholtz equation given by mutually 
orthogonal and different in parity solenoidal vector fields - the poloidal and the toroidal 
fields of fundamental basis (Chandrasekhar, 1961). The fluctuating displacements in 
spheroidal mode are given by the polar (normal or even parity) poloidal vector field 
u s = V x V x \rji(kr)Pi(cos9)] expi(ut) and in torsional mode by the axial (abnormal 
or odd parity) toroidal field u t = V x [r^(A;r)P^(cos 0)\ expi^out). Here je(kr) stands 
for the spherical Bessel function and Pe(cos8) for the Legendre polynomial of degree t. 
Based on this it was shown that eigenfrequencies of spheroidal and toroidal modes can be 
computed from the boundary condition of stress free surface: n k pi k \ r= R = 0, with n k being 
components of the unit vector normal to the globe surface. This condition, applied to each 
of the above two fields of material displacements, leads to the transcendent dispersion 
equations for spherical Bessel functions and its derivatives, F(je(z),dje(z)/dz)\ r= R = 
where z = kr. The roots zu = k^R of these equations uniquely define the frequency 
con = c t ku = z ei c t /R where the subscript i stands for the node number in overtone of 
given multipole degree £. 
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The oscillatory behavior of a perfectly elastic solid sphere in the regime of standing 
waves has been the objective of extensive study in different areas of basic and applied 
researc and perhaps most thoroughly in geophysical context of free oscillations of the 
Earth (e.g. Lamb 1945; Jeffreys 1970; Lapwood and Usami 1981; Aki and Richards 
2002). In this paper the emphasis is laid on less investigated case of the solid elastically 
deformable globe responding to applied compression free load by local fluctuations in 
shear viscoelastic stresses obeying the conditions VkPik = and V k 7i ik = 0. This spe- 
cific case is of particular interest for several problems of nuclear collective dynamics and 
asteroseismology of neutron stars. 



3 Quasistatic oscillations of a viscoelastic solid globe 

The characteristic feature of the solid globe response by small amplitude fluctuations in 
shear elastic p ik = 2/x ^ and viscous 7i ik = 2r] u ik ^ stresses obeying the conditions, 
VfcPifc = and V&TTjfc = 0, is that the solenoidal field of material displacements is subjected 
to the vector Laplace equation 

V 2 u(r,t) = V-u(r,t) = 0. (12) 

In the above context, this last equation can be thought of as long- wavelength limit of the 
vector Helmholtz equation for the standing-wave regime: in the limit of extremely long 
wavelengths A — > oo the wave vector k = (2tt/X) — > 0. 

In order to see what goes on with the energy in the process of such fluctuations 
of material stress, consider the integral equation of energy balance of viscoelastic shear 
deformations 

dV = - p ik ii ik dV- / 7T ik u ik dV (13) 



-2 / (j, u ik ii ik dV - 2 / rju ik u ik dV (14) 



dt J 2 



which is obtained on integrating ([7j) over the globe volume and taking into account the 
condition of stress free surface implying fi\ r= R = and r]\ r= R = 0. One sees that the 
integrands in the right hand side of (Tl4"|) . defining the work done by elastic and viscous 
stresses are independent of spatial derivatives of neither viscoelastic parameters /i and 
v nor shear strain Ui k and rate-of-strain iii k . These latter tensors enter these integrals 



3 It may be relevant to note that in current investigations of technological interest, the model of vibrat- 
ing solid sphere is invoked to explain spectra of ultrafine nanoparticles of ordinary, non-quantum, solids 
(e.g. Hernandez-Roses, Picquart, Haro-Poniatowski, Kanehisa, Jouanne and Morhange 2003; Duval, 
Saviot, Mermet and Murray, 2005; Bastrukov and Lai 2005). 
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as quadratic combinations of the form Ui k Uik and Ui k Ui k , that is, the indegrands do not 
comprise second order spatial derivatives of Ui obeying the Laplace equation (fl2l . In the 
meantime, the bulk forces of both the Hookean elastic stress //* and the Newtonian viscous 
stresses turn to zero: fP = V k Pik = /iV 2 Wj — > and = Vk^ik = rjS/ 2 Ui — > 0. From 
this the term quasistatic, that is force-free, shear oscillations is derived. This observation 
suggests that analysis of quasistatic oscillations can be performed by working from the 
equation of energy balance, not form dynamical equation as is the above outlined case of 
an elastic solid globe responding to perturbation by standing waves of elastic shear. In 
doing this we take advantage of the energy variational method which is due to Rayleigh 
(e.g. Lamb 1945; Lapwood and Usami 1981). 

3.1 The energy variational method 

The point of departure is to use separable representation of the vector field of displace- 
ments and tensor field of shear deformations 

Uj(r, t) = Oi(r) oc(t) u ik (r, t) = a ik (r) a(t) a ik =^(V i a k + V k ai) (15) 

where a(t) is the temporal amplitude and a(r) is the solenoidal vector field of instanta- 
neous, time- independent, displacements. On substituting (fT5|) in (TT4"]) we arrive at the 
well-familiar equations for the amplitude a(t): 

OH r Ma 2 Ka 2 _ Da 2 , , 

- = . 2T H= — + — T = — (16) 

Ma + Da + Ka = (17) 

describing damped harmonic oscillator. The Hamiltonian 7i stands for the total energy 
of dissipative free, normal, vibrations and T is the Rayleigh's dissipative function. The 
integral coefficients of inertia M, stiffness K and viscous friction D are given by 



M = J p(r) ai a t dV K = 2 J fi(r) a ik a ik dV D = 2 J r)(r) a ik a ik dV. (18) 
The solution of (IT7|) taken in the form a(t) = a exp(At) leads to 



A = — — ±i 
2M 



K D 2 



M AM 2 



1/2 

= -r- 1 ±iQ (19) 



where r is the lifetime and Q is the frequency of oscillations damped by viscosity. As a 
result we obtain a(t) = «o exp(— t/r) cos(fit) with 

Q 2 = u 2 [l - (ur)- 2 ] u 2 = ^ r= 2 -^. (20) 
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By 10 is denoted the frequency of the free, non-damped, normal oscillations whose timing 
is determined by the condition ujt » 1. Thus, to compute the frequency and lifetime of 
quasistatic oscillations one need to specify the fields of instantaneous material displace- 
ments a(r) entering the integral coefficients M, K and D of oscillating solid globe. The 
radial profiles of density p(r), the shear modulus /i(r) and shear viscosity r](r) in the solid 
globe are regarded as input data of the method. 

3.2 Material displacements in spheroidal and torsional modes 
of quasistatic oscillations 

It is clear from above that the fields of instantaneous material displacements a(r) obey 
too the vector Laplace equation 

V 2 a(r) = V-a(r)(r) = 0. (21) 

Adhering to the above Lamb's classification of the vibrational eigenmodes of a perfectly 
elastic solid sphere, as spheroidal (shake or s — mode) and torsional (twist or t — mode), 
in the case under consideration, the eigenmodes of quasistatic regime of oscillations can 
be specified by two fundamental solutions to the vector Laplace equation (I2ip . built on 
the general, regular in origin, solution of the scalar Laplace equation 

V 2 x(r) = x(r) = r l P e (cos 9). (22) 

The first one, describing instantaneous displacements in spheroidal mode of quasistatic 
oscillations is given by the even parity (polar) poloidal vector field 

a.= ^VxVx[r X (r)]. (23) 

Remarkably, 

V x V x [r/P*(cos0)] = (£+ l)V/P £ (cos#) (24) 

and, hence, f[2"B"j) can be represented in the following equivalent form 

a s = A s Vr e P e (cos6) (25) 

exhibiting irrotational character of quasistatic spheroidal oscillations: V x a s = 0. For 
the further purpose we note that the velocity field of material displacements 

<5v s (r, t) = u s (r, t) = a s (r)a(t) (26) 
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in the dipole overtone of spheroidal mode is uniform and thereby corresponds to the 
center-of-mass motion with no change of intrinsic material stresses both elastic pij. and 
viscous 7Tjfc. This means that poloidal dipole field, a s (£ = 1) = constant, describes motions 
of translation character, not oscillatory. 

The second fundamental solution to the vector Laplace equation describing instan- 
taneous displacements in the torsional mode is given by the odd parity (axial) toroidal 
vector field 

a t = A t V x [r X (r)] = A t [V X (r) x r]. (27) 

The velocity of material displacements in the toroidal mode Sv t (r, t) = Ut(r, t) = a t (r)d(t) 
can be represented in the well-familiar rotational form 

<fv t (r, t) = [6u(t, t) x r] 5u(r, t) = x 5v t (r, t)} (28) 

exhibiting differentially - rotational, vortical, character of torsional oscillations with the 
angular velocity 

5w(r,t) = A t V r l P e (cosO)a(t). (29) 

This representation of toroidal solution to the vector Laplace equation expresses natural 
extension of the kinematics of rigid-body or uniform rotation to the case of differential 
rotation: the angular velocity is the vector field, not a constant vector. The case of 
rigid-body rotation about polar axis z 

5v = [5u) x r] 5uj = -V x 5v = constant (30) 

is the particular dipole (£ = 1) case of toroidal field given by fTSBl and fT29l . The radial 
profiles of these fields have no nodes in volume of oscillating globe, that is, in the interval 
[0 < r < R]. In view of this, the term non-radial shear oscillations to this type of 
oscillatory behavior is applied. The intrinsic distortions in a viscoelastic sphere undergoing 
non-radial spheroidal and torsional quasistatic oscillations are shown Fig.l. 

One of the important features of the energy variational method is that it can be used 
for a wide class of problems of different physical interest dealing with non-radial vibrations 
of spherical masses of viscoelastic matter characterized by non-homogeneous distribution 
of bulk density and non-uniform profiles of shear modulus and shear viscosity, provided 
that the fields of induced material displacements in the oscillating object obey the vector 
Laplace equation. 
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Spheroidal mode 
L=2 L=3 




Toroidal mode 



L=2 L=3 




Figure 1: Artist view of intrinsic distortions in a viscoelastic solid globe undergoing qua- 
sistatic shear oscillations in quadrupole (L=2) and octupole (L=3) overtones of spheroidal 
and toroidal modes. 
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3.3 Homogeneous solid globe model of uniform density and con- 
stant transport coefficients of viscoelasticity 

For the solid globe with p = const, p = const, rj = const the parameters of inertia M, 
stiffness K and friction D in spheroidal mode of quasistatic oscillations are given by 

M s = AitpA 2 s R 2l+l 1 ^- l K s = &irpA 2 s R 2l - 1 i{l-l) (31) 
D s = 87rr ] A 2 s R 2e - 1 e{£-l). (32) 

One sees, the arbitrary constant A s enter the last equations as A 2 and, hence, the fre- 
quency and lifetime of s — mode, which are determined by ratios of these parameters, 
are independent of specific form of arbitrary constant A s ^ 0. Equations fl3Tl) and (1321) 
show that monopole £ = vibrations cannot be excited due to adopted approximation 
of non-compressional deformations. The dipole field of material displacements describes 
translational motion of the globe as a whole with constant kinetic energy. In this case the 
integral coefficients of both the elastic restoring force, K s , and the damping friction force, 
D s vanish. Thus, the quadrupole I = 2 overtone of spheroidal of quasistatic oscillations 
is the lowest one. The spectral formulas for the frequency of dissipative free oscillations 
uj 2 = K s /M s and the time of their viscous damping r s = 2M S /D S as functions of multipole 
degree i are given by 

a* = a* P(« +!)(/-!)] r. = ^ - ^ - 1} (33) 

2 V- P R2 fo A \ 

W ° = ~p& T ° = — (34) 

where ujq is the natural unit of frequency and Tq of the lifetime of shear vibrations. 

Computed in a similar fashion the integral parameters of torsional non-radial oscilla- 
tions read 

M t = AnpR 2£+ \ /^f/^ ^ K t = An fl R 2t+1 ^f~ 1) (35) 
H (2£+l)(2£ + 3) P 2£+l K ' 

D t = 47rr ] R 2e+l A2A f' 1) . (36) 

The eigenfrequency of non-dissipative oscillations u 2 = K t /M t and the time of their 
viscous dissipation r t = 2M t /D t are given by 

= + 8)(*-l)] Tt = w ^_ ) (37 ) 

where c^o and To are the above defined natural units of frequency and lifetime of these 
fundamental modes. From spectral formulas (|33|) and (|37|) it follows: the larger multipole 
degree of vibration I the higher frequency and the less lifetime, as pictured in Fig. 2. 
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Figure 2: Frequency and time of viscous damping as functions of multipole degree L for 
quasistatic spheroidal and torsional non-radial oscillations of a homogeneous viscoelastic 
solid globe in natural units of the frequency of uj 1 = /i/ (pR 2 ) and the lifetime r = pR 2 /i] 
of shear vibrations. 
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Its worth noting that taking A t = R~( e ~^ we find that at £ = 1 the vector of vorticity, 
5u> = V x 5v = e z d, is uniform vector directed along polar axis. With this A t , the above 
computed parameter of inertia M t as a function of multipole degree £ becomes 

M >M = ^ R * WTwh - (38 > 

It is easy to see that at £ = 1 this parameter equals to the moment of inertia of rigid 
sphere of mass M.: 

2 An 
M t {£= 1) = J= -MR 2 M = — P R 3 . 

5 3 

This is the case when solid globe sets in rigid-body rotation with constant kinetic energy. 
In this case the spring constant K t , and parameter of viscous friction, D t , turn to zero 
and the total energy equals the kinetic energy of rigid body rotation: T = Juo 2 /2. Thus, 
the torsional quadrupole £ = 2 vibration, the twist, is the lowest overtone of the toroidal 
oscillatory mode too. 

It may be noted that frequency of quasistatic oscillations of a perfectly elastic solid 
globe can be derived by staring from quite different footing (Bastrukov 1994), and as 
was shown by Zhen Ye (2000) these spectral formulas can be regained by working from 
solutions of the standing wave regime of oscillations by taking in resultant dispersion 
equations, for both spheroidal and toroidal modes, the long wavelength limit. However, 
this latter approach becomes less efficient when shear modulus is a function of position. 
Also, it is notable that the effect of viscosity on the motion of solid material continuum 
is described by the same Newtonian law as in the viscous liquid (Landau, Lifshitz, Ko- 
sevich and Pitaevskii 1995). Thereby it seems worthy of comparing the above results for 
the damping time of quasistatically oscillating viscoelastic solid globe with that of fluid- 
mechanical calculation for time of viscous damping of non-radial oscillations of a liquid 
drop. 



3.4 Fluid-mechanical approach to the viscous damping of non- 
radial oscillations of an incompressible liquid drop 

The fluid mechanical problem of damping of irrotational oscillation of a spherical mass of 
an incompressible viscous liquid has first been tackled by Lamb (1945) too based on the 
Navier-Stokes equation of hydrodynamics of viscous incompressible liquid 

p5vi = -Vi5p + V k 5n ik 5n ik = 2r]5v ik (39) 
Sv ik = -[Vi5v k + V k 5vi\ Sv kk = V/c Sv k = 0. (40) 
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The equivalent vector form of (J3HD reads 

p 5v = -V5p + r]V 2 5v V • 5v = V x 5v = 0. (41) 

Note that irrotational character of the field of velocity <5v in fluctuating incompressible flow 
implies 5vi = — V* Sx, and hence: V • 5v = — > V 2 6x = where 5\ = Ai r Pt(cos 6) at(t), 
so that <5v is the solution to the vector Laplace equation V 2 5v = 0. Understandably 
that the condition of irrotational flow (V x 5v = 0), typical of hydrodynamical consider- 
ation, excludes solution of the vector Laplace equation describing differential rotational 
fluctuations. In Lamb's calculations by Rayleigh' method the focus was laid on the non- 
compressional oscillations in which 5p = 0. In this case the energy balance equation 



di I ^T dV = ~ I 67lik6vikdV = ~ 2 V^ ! ikSv ik dV. (42) 



with help of substitution 5vi(r, t) = a«(r) a(t) is reduced to 

— = -IT T = T = 43 

dt 2 2 v ; 

M = Jpa^dV D = 2 J 7] a ik a ik dV a ik = -(Vja fe + V fc aj). (44) 

where T is kinetic energy of oscillations and JF is the dissipative function. Note, the 
analytic form of inertia M and viscous friction D is precisely coincides with the above 
obtained integral parameters of quasistatic oscillations of a viscoelastic solid globe. The 
next step is to represent temporal amplitude of oscillations in the form 

a(t) = A(t) cos ujt (45) 

where uj is frequency of oscillations and A(t) smoothly varying function of time. With 
this form of a the mean value (averaged over the period of vibration) of the kinetic 
energy and dissipative function are given by < T >= (l/2)u 2 M A 2 (t) and < T >= 
(1 / 4) D u 2 A 2 (t) , respectively. The key assumption of this line of argument is that the 
equation of energy balance, (|43p . expressing the link between the time-dependent kinetic 
energy T and doubled dissipative function 2J 7 , preserves its validity for the mean values 
of these quantities 

^<T>=-2<T> - M^- = --DA(t) (46) 
dt dt 2 w y J 

2M 

A(t)=A exp(-t/r) r = — . (47) 

Thus, the computational formula for the viscous damping time r is identical to that ob- 
tained for a viscoelastic solid globe. This explains why the spectral formula for the time 
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of viscous damping of irrotational quasistatic spheroidal oscillations of the solid globe 
exactly coincides with the Lamb's formula for the damping time of the liquid drop oscilla- 
tions. In the context of viscous damping of global non-radial pulsations of neutron stars, 
the above spectral formulas for lifetime of spheroidal and toroidal modes are discussed 
in (Cutler and Lindblom 1987; Bastrukov, Weber and Podgainy 1999; Bastrukov, Yang, 
Podgainy and Weber 2003). 

3.5 Nuclear giant resonances in the model of quasistatic oscil- 
lations of viscoelastic solid globe 

The physical idea underlying the nuclear solid globe model is that the order encoded 
in the shell-model treatment of nuclear structure signifies quantum solidification of de- 
generate Fermi-system of nucleons. In this model a nucleus is thought of as an ultrafine 
(femtometer dimension) particle of Fermi-solid - the quantum (fermionic) condensed mat- 
ter endowed with material properties of viscoelasticity typical of ordinary solids. As was 
mentioned, the credit for recognition of solid mechanical elasticity in the nuclear giant- 
resonant response goes to seminal work of Bertsch (1974). The basic inference of this 
influential paper is that experimentally established smooth dependence of the energy of 
isoscalar giant quadrupole resonance can be understood as manifestation of spheroidal 
mode of non-compressional, that is pure shear, quasistatic oscillations of an ultrasmall 
spherical particle of an elastic Fermi-solid whose shear modulus equals to the pressure of 
Fermi-degeneracy: /x = Pp = pvp/5 (see also, Semenko 1977; Nix and Sierk 1980; Wong 
and Azziz 1981; S. Mi§icu 2006). By now there are quite cogent arguments showing that 
both the giant electric and magnetic resonances can be treated on equal solid-mechanical 
footing as manifestation of spheroidal (electric) and torsional (magnetic) oscillations of 
nuclear femtoparticle (e.g. Bastrukov, Mi§icu and Sushkov 1993). The electromagnetic 
nomenclature of nuclear giant resonances owe its origin to the type of induced electromag- 
netic moment of fluctuating current density: the giant-resonance excitations of electric 
type are associated with spheroidal mode of quasistatic oscillations of irrotational field 
of material displacements u s and those for the magnetic type with torsional mode of 
quasistatic differentially-rotational oscillations of toroidal field of material displacements 
u t . The multipole degree A of electric Ai(EX) and magnetic M.(MX) moments equal to 
multipole degree I of excited spheroidal and torsional oscillatory state, respectiveljfl. 
Fig. 3 shows how the spectral formulas fl33|) and (j3~T|) of the nuclear solid globe 

4 The electric and magnetic multipole moment of electric current density are given by 
M{E£) = N Ee [ j(r, t) • V/ P t dV M(M£) = N Ml f j(r, t) • [r x V] / P t dV 
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Figure 3: The energy E = Huj and width T = hr^ 1 of isoscalar electric and magnetic 
resonances, interpreted from equal footing as manifestation of spheroidal (electric) and 
torsional (magnetic) quasistatic oscillations whose frequency and lifetime is given by ( 1331) 
and ( |37l) . respectively. The multipole degree A of excited electromagnetic moment of 
the electric current density equal to multipole degree of oscillations A = i > 2, and 
ujp = vp/R, as pointed out in the text. 
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model can be utilized to extract numerical estimates for the shear modulus and shear 
viscosity of nuclear matter from the experimental data of nuclear physics. The nuclear 
density p = 2.8 10 14 g cm~ 3 and the nucleus radius R = r A 1 ^ (r = 1.2 1CT 13 cm) are 
well defined quantities. Therefore by fitting the data on general trends in the energy 
centroids and spread widths of isoscalar (with isospin quantum number T = 0) giant 
resonances, which are computed with use of standard quantum mechanical equations 
for the energy centroid E = hu and spread width T = h/r, we can extract the values 
of two unknown quantities - the shear modulus /i and shear viscosity rj of the isospin- 
symmetric nuclear matter, as was mentioned with T = 0. Note, taking into account that 
in this model p = (l/5)pv F , the natural unit of frequency of shear elastic oscillations u 
carrying information on shear modulus can be expressed in terms of the Fermi-frequency 
oop = Vp/R as follows = p/(pR 2 ) = (l/5)up. It is to be emphasized that in the 
presented in Fig. 3 comparison of theoretical calculations (lines) with data (symbols), 
the transport coefficients of solid-mechanical elasticity p and 77 have been fixed so as to 
reproduce overall trends in data on the energy centroid and the spread width of solely 
one type of resonance - the giant electric quadrupole (E2) resonance associated with 
quadrupole overtone of spheroidal quasistatic oscillations. The predictive power of the 
solid globe model is demonstrated by fairly accurate account of experimental systematics 
on the electric octupole (E3) and the magnetic quadrupole (M2) giant resonances which 
is attained with no use of any adjustable constants. 

As was mentioned, the giant magnetic resonances exhibiting torsional oscillations of 
nucleus is the mode of giant-resonance nuclear response that can be understood only 
on the basis of the solid globe model, not a liquid drop. One more demonstration of 
the predictive power of the nuclear Fermi-solid globe model is pictured in Fig.4 showing 
adequate description of the electromagnetic excitation strength (total excitation probabil- 
ity) - B(M2) of giant magnetic resonance versus mass number A which is attained with 
no use of any adjustable constants (Bastrukov, Molodtsova and Shilov 1995). Also, its 

where j(r, t) = p e u(r, t) with p e = (Z/A)n is the charge density and n being the nucleon density. The ex- 
act expressions of Nei and Nmi can be found elsewhere (e.g. Alder and Steffen 1975; Greiner and Maruhn 
1996). Note, the quasistatic oscillations of irrotational field of material displacements in spheroidal mode 
u = u s = a s (r) a(t) result in excitations of vibrational states with non-zero electric multipole moment, 
whereas magnetic multipole moment for this kind of oscillations is zero. The opposite situation takes 
place for quasistatic oscillations of differentially-rotational displacements u = u t = a t (r)a(i) in tor- 
sional mode which lead to vibrational excited states with non-zero magnetic multipole moment, while 
the electric multipole moment is zero. The computational formulas for the total excitation strength 
(probability) of giant electric and magnetic resonance are given by B(E£) = (2£ + 1) < \M(El)\ 2 > and 
B(M£) = (2£+ 1) < \M(Ml)\ 2 >, respectively, where bracket stands for time average. 
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Figure 4: Theoretical calculation (lines) for general trends throughout periodic table of 
the energy centroid and the total excitation probability (computed with no adjustable 
parameters) for the nuclear giant magnetic resonances of multipole degree A > 2 treated 
in the from stand point of the Fermi-solid globe model as manifestation of torsional 
quasistatic oscillations versus experimental data (symbols) of S-DALINAC data (Luttge, 
von Neumann-Cosel, Neumeyer and Richter 1996). 
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worth emphasizing that the spectral formulas for the quasistatic shear oscillations match 
much better the experimental data as compared to calculation implying oscillations in 
the regime of standing waves, extensively studied in (Wong and Azziz 1981). 

The above procedure of extracting the values of transport coefficient of solid-mechanical 
viscoelasticity of nuclear matter - the shear modulus p and the shear viscosity 77 yields: 

p — 4 • 10 33 dyn cm -2 rj = 3 • 10 11 dyn seccm" 2 

(see also Hasse 1978; Nix and Sierk 1980; Kolomietz 1990). For comparison, the density 
of steel p ~ 10 g cm -3 and the shear modulus p ~ 10 11 dyn cm" 2 (e.g. Fetter and 
Walecka 1980). In theoretical geophysics this last estimate is of the same order as the 
average value for the shear modulus of the Earth's material (e.g. Lapwood and Usami 
1981). It is appropriate mention here that obtained in this manner estimate for the 
shear viscosity is much less than that derived in (Flowers, Itoh 1976: rj ~ 10 18 dyn sec 
cm" 2 ) by use of absolutely different theoretical arguments (see also Cutler and Lindblom 
1987; McDermott, Van Horn and Hansen 1988; Bastrukov, Yang, Podgainy and Weber 
2003). This inference of nuclear physics is crucial to the asteroseismology of neutron 
stars studying the timing of oscillatory behavior of pulsars and magnetars. With the 
obtained from giant resonances value of shear viscosity 77, one finds that the time of 
viscous damping of shear oscillations, r = pR 2 /i], in the object of nuclear density p 
and radius R ~ 10 6 cm, typical of neutron stars, is evaluated as r ~ 10 7 — 10 s years. 
The net outcome of this finding is that neutron stars are expected to be in the restless 
state of background shear vibrations (triggered either by supernova event or starquake) 
during all its life as a pulsar. As is shown in the following section this inference is in 
line with astrophysical observations and conclusions of several independent theoretical 
investigations on the oscillatory behavior of neutron stars. 

4 The solid globe models of astrophysical interest 

The above homogeneous models of viscoelastic solid globe and viscous liquid drop are 
of practical usefulness when there are physical arguments justifying approximation of 
uniform profiles for the density and the transport coefficients of shear viscoelasticity. This 
approximation surely becomes unwarranted when studying global vibrations of a highly 
stratified in density solid mass, as is the case of solid Earth-like planets and collapsed 
degenerate stars, like white dwarfs and pulsars. In this case the shear modulus and shear 
viscosity cannot be considered as having constant values throughout the star volume (e.g. 
Jeffreys 1976; Lapwood and Usami 1981). It is the purpose of this section to investigate 
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quasistatic oscillations in several mathematical models of viscoelastic solid globe with 
non-uniform profiles of density and parameters viscoelasticity. 

Bearing in mind the astrophysical context of the models under consideration, in what 
follows we discuss quasistatic global shear vibrations of self-gravitating masses whose 
state of gravitational equilibrium is determined by coupled equations for the potential of 
self-gravity and equation for the pressure 

V 2 $(r) = -47rGp(r) VP(r) = p(r) V$(r). (48) 

which, of course, valid for both liquid and solid stars, by virtue of universal character of 
Newtonian law of gravitation. In these equations, the density profile p(r) is regarded as 
a given function of distance from the star center. Then having solved the first equation 
in this line for the potential of self-gravity $(r) we can obtain the profile of pressure 
P(r) providing gravitational stability of the star. In doing this the standard boundary 
conditions of gravitostatics for the potential on the globe surface, Q in = $ ou t\r=R and 
V r $j n = V r & out \ r =R, and the condition of stress-free surface for the pressure, P r =B, = 0, 
should be used. 

As a representative example, for the standard model of homogeneous star with p = 
constant throughout the star volume, from (j48l it follows that profile of pressure inside 
the star is given by P(r) = P c [l — (r/R) 2 ], where P c = (2n/3)Gp 2 R 2 is the pressure in the 
star center. On the other hand, it should be clearly realized, the central pressure must be 
evaluated by use of physical arguments regarding the equation of state of stellar matter 
and dynamical processes opposing the self gravity. For instance, in a highly idealized 
model of neutron star, thought of as a homogeneous spherical mass of self-gravitating 
degenerate neutron matter condensed to the normal nuclear density, the central pressure 
counterbalancing the self-gravity is taken in the form of equation of state for degenerate 
Fermi-gas of non-relativistic neutrons Pp = Kp 5 ^ 3 . Then, equating the central pressure 
P c = {2ir/3)Gp 2 R 2 obtained from the equation of gravitational equilibrium (|48|) with Pp 
we obtain the equation for the neutron star radius R which leads to canonical estimates 
R 9a 12 km and mass M = (4ir/3)pR 3 ~ 1.3 M . This is perhaps the most simple ex- 
ample showing that the neutron star radius is primarily determined by the behavior of 
degeneracy pressure of the neutron star matter in the vicinity of the nuclear saturation 
density (e.g. Lattimer and Prakash 2001). Also, it may be worthy of noting that from 
quantum-mechanical side, the fact that the pressure in the star counterbalancing New- 
tonian self-gravity is determined by the degeneracy pressure of Fermi-gas of independent 
neutron-like quasiparticles means that the single-particle states of incessant quantum- 
wave Fermi motion of an individual quasiparticle in the mean field of self-gravity are 
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described by Schrodinger equation coupled with Poisson equation for the potential of the 
mean gravitational field $(r) in the star 

if) U(r) = — m n $(r) 

777 

V 2 $(r) = -4%Gp(r) p = — £ k% = 2.8 x 10 14 g cm" 3 (49) 

where U is the potential energy of the in-medium neutron quasiparticle in its unceasing 
Fermi-motion whose collision free character is provided by Pauli exclusion principle; p 
stands for the uniform density in the star bulk equal to the normal nuclear density and 
m n is the effective mass of the neutron quasiparticle. The solution to (|49|) having the form 
$(r) = (27r/3) Gp (3i? 2 — r 2 ) suggests that the potential energy U can be represented in 
the well-familiar form of spherical harmonic oscillator 

U = -m n ^r) = -U G + ^^- U G = 2irmpGR 2 u G = ^-Gp (50) 

where U G is the depth of spherical gravitational trap and uj g is the characteristic unit 
of frequency of gravitational vibrations. The stationary states of single-particle Fermi- 
motion of individual neutron in the mean field of self-gravity are described by Hamiltonian 

if ) = eif) e = E + U G (51) 

whose spectrum of energies, accounted from the bottom of potential well, is well-known 
e^v = £o(N + 3/2) where eo = hue ~ 10~ n eV is the energy of zero-point oscillations which 
is the measure of the energy distance between single-particle states of neutron quasipar- 
ticles in the potential of self-gravity of neutron star (note the average distance between 
discrete states of degenerate electrons in terrestrial solids such as metals and semiconduc- 
tors A ~ 1CT 18 — 10~ 20 eV). Here N = 2£ + n is the shell's quantum number of harmonic 
oscillator, n and i are the principle and orbital quantum numbers of single-particle orbits, 
respectively. This shows that the shell-ordered clusterization of single-particle energies 
of neutron quasiparticles in the potential of mean gravitational field of the neutron star 
model has common features with that for nucleons in the mean field potential of the 
nuclear shell model having different physical origin. Thus, in addition to the similar- 
ity in nucleon-dominated composition and the fact of common genetic origin of nuclei 
and neutron stars rooted in supernovae of second type, an identical quantum-mechanical 
mean field treatment of orderly-organized intrinsic single-particle structure of these ob- 
jects indicates that they are made of the identical in its material properties degenerate 
Fermi-matter endowed with solid-mechanical property of viscoelasticity - elastically de- 
formable Fermi-solid. 
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One more cognitive value of this historically first and admittedly idealized homoge- 
neous model (using the only one input parameter - the density of nuclear matter) is in 
highlighting the compressional effect of gravity, that is, in showing that material pressure 
is increased in the direction to the star center whereas the density of stellar matter re- 
mains constant (the well-familiar example is pressure increase with the depth of ocean 
whereas the density of water is practically unaffected). 

Bearing this in mind, it seems not inconsistent to assume that compressional effect of 
gravity on the viscoelastic parameters of self-gravitating solid mass is analogous to that for 
the pressure P. More specifically, it is assumed that the shear modulus profile p, = p(r) 
(having physical dimension of pressure P) and the shear modulus profile 77 = r](r) are 
identical to that for the pressure profile P(r). Based on this conjecture, in the reminder 
of this paper we compute the frequency and lifetime of quasistatic oscillations for several 
models of astrophysical interest by evaluating first the radial profile of pressure from 
equations of gravitational equilibrium (|4"5j) . 

4.1 Solid star model with uniform density and non-uniform pro- 
files of shear modulus and shear viscosity 

Following the above line of argument, consider first the model of uniform density and 
non-uniform shear modulus and shear viscosity whose profiles are taken to be identical to 
that for the pressure of self-gravity: 



where p c and u c the shear modulus and shear viscosity of stellar matter in the star center, 
respectively. Then, for integral parameters of global quasistatic oscillations in spheroidal 
mode we obtain 



p = constant 




(52) 




(53) 



(54) 



and for the frequency and lifetime we get 




pR 2 



(55) 
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The integral parameters of quasistatic differentially rotational oscillations are given by 



M '-^ A ' g2M wrmh (56) 
* = ^ RWA 'iSmrT) (57) 

* = ***^* { u?iM + s> (58) 

and for the frequency and lifetime of toroidal mode we obtain 

^ = 2^(1-1) r t = j^. (59) 

Thus, the compressional effect of gravity on transport coefficients of solid-mechanical 
viscoelasticity of stellar matter is manifested in the spectral formulas for the frequency 
and lifetime of global oscillations. Specifically, for p c and r\ c taken equal to p and 77 of the 
homogeneous model, we find that the frequency of quasistatic global oscillations of the 
non-homogeneous solid globe is lower than that for homogeneous one and, hence, period of 
oscillations, P = 2tt/uj, is longer. The lifetime of such oscillations in a non- homogeneous 
solid globe are longer as compared to that in the model with homogeneous density and 
uniform profiles of parameters of viscoelasticity. 

Fig. 5 shows that periods of background elastic pulsations of a homogeneous neutron 
star model computed with use of spectral formulas (|55p and (|59p with the value of shear 
modulus extracted from data on nuclear giant resonances coincides with the timing of 
microspikes in the windows between the main pulses. Taking into account the above 
inference regarding the damping time of global quasistatic elastic vibrations of pulsars by 
viscosity of neutron star matter one can conclude that these microspikes owe its origin to 
non-radial elastodynamic (EDM) spheroidal and torsional pulsations triggered by neutron 
star quakes (e.g. Van Horn 1980, McDermott, Van Horn and Hansen 1988; Bastrukov 
Weber and Podgainy 1999). 

Deserved for special comment is the following particular case. As was mentioned, the 
physical dimension of shear modulus p is the same as of pressure P. Bearing this in mind, 
consider the model of homogeneous solid star putting 



p = constant p(r) = p c 



1 --b 



p c = P c = (2n/3)Gp 2 R\ (60) 



In terms of this last parametrization of p c , the eigenfrequencies of global non-radial 
spheroidal and torsional quasistatic oscillations are represented as follows (Bastrukov 
1993, 1996) 

4-7T 

ul = 2ul{l-l) J* = <4(l - 1) u% = - P G. (61) 
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Figure 5: Periods P against multipole degree L for spheroidal and torsional elastody- 
namic (EDM) modes of background non-radial pulsations of neutron star manifested by 
microspikes of millisecond duration in the windows of the rotation driven main pulse train. 
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Figure 6: Fractional frequency uj/ujg versus multipole degree L of non-radial oscillations 
for fundamental spheroidal and toroidal gravity modes in a solid star in juxtaposition 
with that for the Kelvin's fundamental gravity mode in a liquid star. 



For a solid star - self-gravitating mass of a viscoelastic solid, the last spectral formulas 
seems to have the same physical meaning as Kelvin's formula does 
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for a liquid star - self-gravitating mass of an incompressible inviscid liquid. This last spec- 
tral formula can be obtained by the above expanded method from equations of canonical 
hydrodynamical model for a self-gravitating inviscid liquid relying on equations of fluid 
mechanics coupled with the equation of self-gravity 
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(e.g. Lamb 1945; Chandrasekhar 1961; Bastrukov 1996). The comparison of flBTT) and 
(l62~l) is depicted in Fig. 6. 
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4.2 Solid star model with the non-uniform density and non- 
uniform profiles of shear modulus and shear viscosity 



Consider a solid star model with the non-uniform distribution of density 

2 



p(r) = p c 



1-1- 
R 



(64) 



After long but simple calculation of pressure from equations of gravitational equilibrium 
18]) we obtain 

2 



P(r) = Pr 
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Taking into account the above conjecture of similarity in form of radial profiles for vis- 
coelastic parameters and pressure we put 
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where p c and v c are the sear modulus and shear viscosity in the globe center. Omitting 
tedious calculation of integrals for M, K and D, which has been performed with help of 
results presented in appendices, we obtain 
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To (£-l)(2£+ll) Tt = T °(£-l)(2£+13) T ° = f " (68) 
Note, putting p c of this non-homogeneous model equal to the density of homogeneous 
model of previous subsection p we find that these stellar models have different radius and 
mass. The considered in this subsection a non-homogeneous model of a solid star is bigger 
in radius and more massive in mass, this difference, as is clearly seen, is traced in the 
values of periods and lifetimes for both s-mode and t-mode of quasistatic oscillations. The 
presented calculation justifies inferences of foregoing subsection regarding the frequencies 
and lifetimes of quasistatic spheroidal and torsional non-radial oscillations computed in 
homogeneous and non-homogeneous solid globe models. 



4.3 Self-gravitating mass with non- uniform singular density and 
non-uniform profiles of shear modulus and shear viscosity 

Finally, consider self-gravitating spherical mass of radius R with singular in origin density 
of the form 
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R 
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Unlike the above case of models with regular density profiles, in the last formula the 
constant p can be fixed by condition on the globe surface: p{r = R) = p s = constant, 
not in the center. This kind of mass distribution are currently utilized in astrophysical 
models aimed at the search for dark matter candidates (e.g. de Boer, Sander, Zhukov, 
Gladyshev and Kazakov 2005). This model is interesting in that the total mass 
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has the form identical to that for a spherical mass of a homogeneous density. Following 
the above prescription and computing first the pressure profile inside this object from 
equations of gravitational equilibrium, we obtain 
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It is seen that the pressure is linear function of distance from the center and in the center 
the pressure is finite. Keeping in mind our conjecture that profiles of shear modulus and 
shear viscosity follow the same radial dependence as the pressure profile does we put 
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For integral parameters of quasistatic oscillations in spheroidal mode we obtain 
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The frequency and lifetime are given by 
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For torsional quasistatic vibrations one has 
M t = A7cpA 2 R 2£+3 - 
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For the frequency and lifetime of toroidal mode we get 
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All the above models preserve one important feature of homogeneous model regarding the 
quadrupole (£ = 2) degree of lowest overtone of spheroidal and torsional quasistatic shear 
oscillations. 
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5 Summary 

A fundamental understanding of oscillatory modes in a spherical mass of viscoelastic solid 
is needed in different areas of physics and astrophysics embracing wide spatial scale of 
material objects ranging from ultrafine pieces of solid matter such as metallic nanopar- 
ticles and atomic nuclei to huge astrophysical objects like planets of the solar systems 
and collapsed degenerate compact stars like white dwarfs and pulsars. In this paper, 
the regime of quasistatic shear vibrations of viscoelastic solid globe has been investigated 
in some details whose non-trivial feature is that the restoring force of elastic stress and 
dissipative force of viscous stress entering the equation of solid-mechanics turn to zero, 
whereas the work done by the above stresses in the bulk of oscillating solid globe does 
not. It is shown that in this regime the frequency and lifetime of quasistatic spheroidal 
and torsional vibrations can be computed working from energy balance equation under- 
lying the Rayleigh's energy variational method. The efficiency of this method has been 
demonstrated by analytic results for several astrophysical models illustrating the effect 
of self-gravity on profiles of solid-mechanical transport coefficients of shear viscosity and 
elasticity and its manifestation in the spectra of quasistatic non-radial oscillations. The 
practical usefulness of considered mathematical models and obtained analytic estimates 
is that they are quite general and can be applied to different spherical systems whose 
behavior is supposedly governed by equations by solid-mechanics. 

6 Appendices 

For computational purpose it is convenient to represent the components of tensor of shear 
deformation 
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in spherical polar coordinates making use of angle variable ( = cos 6. The explicit form 
of these components reads 
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For the spheroidal mode of quasistatic shear oscillations the spherical polar components 
of poloidal field of instantaneous displacements are give by 
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The components of the shear deformation tensor in spheroidal mode of quasistatic oscil- 
lations are given by 
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After integration over the solid angle for parameters of inertia, shear stiffness (and shear 
friction) we obtain 

M s = J p{r)oi a4V = J p{r)[a 2 r + a 2 g ]dV = Aix£A 2 s J pr 2e dr 
K s = 2 J n{r)a ijaij dV = 2 J /i(r) (a 2 rr + a 2 ee + a 2 ^ + 2a 2 e ) dV 
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and result for D s is similar to that for K. The finale computational formulas for frequency 
and lifetime of spheroidal fundamental mode of quasistatic shear oscillations reads 
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In the toroidal fundamental mode of quasistatic shear oscillations, the spherical polar 
components of toroidal instantaneous displacements have the form 
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The non-zero components of strain tensor in torsional mode are given by 
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The inertia, stiffness and friction in torsional mode of differentially rotational non-radial 
shear oscillations after integration over the solid angle have the form 

M t = Jp(r)a t a t dV = J p^dV = AnA 2 ^^ £p(r)r 2e+2 dr 
K t = 2 J p{r) aijUij dV = 2 J p(r) (a^ + a 2 9(f) ) dV 
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The integral coefficient of friction D in torsional mode is computed in a similar fashion. 
Finale computational formula for the frequency and lifetime of toroidal fundamental mode 
of quasistatic shear oscillations is given by 
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The practical usefulness of the calculus presented in this appendices is that it can be 
utilized in the study of diverse class of solid globe models of physical interest. 



Acknowledgment 

This work is partially supported by NSC of Taiwan, Republic of China, under grants 
NSC 96-2811-M-007-001 and NSC-95-211-M-007-050 and under protocol JINR (Russia)- 
IFIN-HH (Romania), project number 3753-2007. 



References 

[1] H. Bethe, Nuclear Physics, Rev. Mod. Phys. 71 (1999) S6-S15. 

[2] R.A. Guyer, The Fermi solid, in Selected Topics in Physics, Astrophysics, and Bio- 
physics (Reidel, Dordrecht-Holland, 1973) pp. 44-49. 



33 



[3] E.Polturak and N.Gov, Inside a quantum solid, Contemporary Physics 44 (2003) 
145-151. 

[4] M. Harwit, Astrophysical Concepts 3th edn (Springer- Verlag, 2006). 

[5] S.L. Shapiro and S.A. Teukolsky, Black Holes, White Dwarfs and Neutron Stars 
(Wiley, New York, 1983). 

[6] F. Weber, Pulsars as astrophysical laboratory for nuclear and particle physics (IOP 
Publishing 1999). 

[7] N. Glendenning, Compact Stars (Springer, Berlin, 2000). 

[8] J.M. Lattimer and M. Prakash, Neutron Star Structure and the Equation of State, 
Astrophys. J. 550 (2001) 426-442. 

[9] F. Weber, Strange quark matter and compact stars, Prog. Part. Nucl. Phys. 54 
(2005) 193—288. 

[10] R.X. Xu, Solid quark stars, Astrophys. J. 596 (2003) L59-L62. 

[11] B.J. Owen, Maximum elastic deformations of compact stars with exotic equations of 
etate Phys. Rev. Lett. 95 (2005) 211101. 

[12] H.M. van Horn, Micropulses, drifting subpulses, and nonradial oscillations of neutron 
stars, Astrophys. J. 236 (1980) 899-903. 

[13] P.N. McDermott, H.M. van Horn and C.J. Hansen, Non-radial oscillations of neutron 
stars, Astrophys. J., 325 (1988) 725-748. 

[14] T.E. Strohmayer, J.M. Cordes and H.M. van Horn, Determining the coherence of 
micropulses, Astrophys. «/., 389 (1992) 685-694. 

[15] S.I. Bastrukov, F. Weber and D.V. Podgainy, On stability of global non-radial pul- 
sations of neutron stars J. Phys. G25 (1999) 107-127. 

[16] D. Pines, J. Shaham and M.A. Ruderman, Neutron Star Structure from Pulsar Obser- 
vations, in Physics of Dense Matter Proc. of IAU Symposium No. 53 (1974) 189-207. 

[17] V. Canuto and S.M. Chitre, Crystallization of dense neutron matter, Phys. Rev. D9 
(1974) 1587-1613 . 

[18] O. Blaes, R. Blandford, P. Goldreich, P. Madau, Neutron starquake models for 
gamma-ray bursts, Astrophys. J. 343 (1989) 839-848. 



34 



[19] L.M. Franco, B. Link and R.I. Epstein, Quaking neutron stars, Astrophys. J. 543 
(2000) 987-994. 

[20] B. Cheng, R.I. Epstein, R.A. Guyer and A. C.Young, Earthquake-like behaviour of 
soft gamma-ray repeaters, Nature 382 (1996) 518-520. 

[21] M. Ruderman, Crystallization and torsional oscillations of superdense stars, Nature 
218 (1968) 1128-1129. 

[22] C.J. Hansen and D.F. Cioffi, Torsional oscillations in neutron star crusts, Astophys. 
J. 238 (1980) 740-742. 

[23] S. Bastrukov, D. Podgainy, J. Yang and F. Weber, Magneto Torsional Pulsations of 
Magnetars, Mem. Soc. Astron. Ital. 73 (2002) p. 522-528. 

[24] Lupin C.-C. Lin and Hsiang-Kuang Chang, Periodicity Search in the X-ray Data of 
RX J0007.0+7302, Astrophys. and Space Science 267 (2005) 361-367. 

[25] T.E. Strohmayer and A.L. Watts, The 2004 hyperflare from SGR 1806-20: Further 
evidence for global torsional vibrations, Astrophys. J., 653 (2006) 593-601. 

[26] T. Strohmayer and L. Bildsten, New views of thermonuclear bursts, in Compact 
stellar X-ray sources Cambridge Astrophys. Series No. 39 (2006) 113-156. 

[27] G.S. Bisnovatyi-Kogan, Dynamic confinement of jets by magnetotorsional oscilla- 
tions, Monthly Notes Roy. Astronom. Soc, 376 (2007) 457-464. 

[28] G.F. Bertsch, Elasticity in the response of nuclei, Ann. Phys. N. Y. 86 (1974) 138-146. 

[29] M.N. Harakeh and A. van der Woude, Giant Resonances: Fundamental High- 
Frequency Modes of Nuclear Excitations (Oxford Science Publishing, Oxford, 2003). 

[30] A. Richter, Giant resonances - experiments at highest resolution, wavelets and scales, 
Prog. Part. Nucl. Phys. 55 (2005) 387-396. 

[31] S.F. Semenko, Collective coordinates in the random phase approximation and the 
systematics of giant multipole resonances, Sov. J. Nucl. Phys. 26 (1977) 143-146. 

[32] H. Sagawa and G. Holzwarth, Giant Resonances: A Comparison between TDHF and 
Fluid Dynamics in Small Amplitude Vibrations of Spherical Nuclei, Prog. Theor. 
Phys. A325 (1978) 1213-1229. 



35 



[33] J.R. Nix and A.J. Sierk, Macroscopic description of isoscalar giant multipole reso- 
nances, Phys. Rev. C21 (1980) 396-404. 

[34] C.-Y. Wong and N. Aziz, Nuclear giant resonances as elastic vibrations, Phys. Rev. 
C24 (1981) 2290-2310. 

[35] R. Hasse, G. Ghosh, J. Winter and A. Lumbroso, Isoscalar giant-resonance energies 
and long-mean-free path nuclear fluid dynamics, Phys. Rev. C25 (1981) 2771-2779. 

[36] S. Stringari, Fluid-dynamical description of nuclear collective excitations, Ann. Phys. 
(N.Y.) 151 (1983) 35-70. 

[37] W. Norenberg, Elastoplasticity of finite fermi systems, Nucl. Phys. 428 (1984) 177- 
187. 

[38] E.B. Balbutsev and I.N. Mikhailov, Dynamics of nuclear integral characteristics, J. 
Phys. G14 (1988) 545-567. 

[39] V.M. Kolomietz, Local density approach to atomic and nuclear physics (Naukova 
Dumka, Kiev, 1990). 

[40] M. Di Toro, Semiclassical approach to the description of collective nuclear motions, 
Sov. J. Part. Nucl. 22 (1991) 185-210. 

[41] S.I. Bastrukov, §. Misicu and A.V. Sushkov, The dipole torus mode in nuclear fluid- 
dynamics, Nucl. Phys. A562 (1993) 191-204. 

[42] §. Mi§icu, Interplay of compressional and vortical nuclear currents in overtones of 
the isoscalar giant dipole resonance, Phys. Rev. C73 (2006) 024301. 

[43] S. Raman, LW. Fagg and R. Hicks, Giant Magnetic Resonances in: Electric and 
Magnetic Giant Resonances ed. by J. Speth (World Scientific, Singapore, 1991) p. 355. 

[44] C. Luttge, P. von Neumann-Cosel, F. Neumeyer and A. Richter, Magnetic dipole 
and quadrupole response of nuclei, supernova physics and in-medium vector meson 
scaling, Nucl. Phys. A606 (1996) 183-200. 

[45] G. Holzwarth and G. Eckart, The nuclear "Twist", Z. Phys. A283 (1977) 219-220. 

[46] S.I. Bastrukov and V.V. Gudkov, The twist MA,T = giant modes in spherical 
nuclei, Z. Phys. A341 (1992) 395-399. 



36 



[47] S.I., Bastrukov, I.V. Molodtsova and V.M. Shilov, Torsional multipole magnetic re- 
sponse of a heavy spherical nucleus, Phys. Scripta 51 (1995) 54-59. 

[48] S.I. Bastrukov, J. Libert and I. Molodtsova, Elastodynamic features of nuclear matter 
from macroscopic model of giant magnetic resonances, Int. J. Mod. Phys. E6 (1997) 
89-110. 

[49] W. Greiner and J. A. Maruhn, Nuclear Models (Springer, Berlin, 1996). 

[50] D.J. Rowe, Nuclear collective motion: models and theory (Methuen, London, 1970). 

[51] H. Lamb, On the vibrations of an elastic sphere, Proc. London. Math. Soc. 13 (1881- 
82) 189-212. 

[52] A.E.H. Love, Treatise on Mathematical theory of elasticity (Dover, New York, 1944). 

[53] R.M. Morse and H. Feshbach, Methods of Mathematical Physics (MCGraw-Hill, New 
York, 1953). 

[54] A.C. Eringen and E.S. Suhubi, Elastodynamics vol. 2 (Pergamon Press, New York, 
1975). 

[55] L.D. Landau, E.M. Lifshits, A.M. Kosevich and L.P. Pitaevskii, Theory of Elasticity 
(Butterworth-Heinenann, Oxford 1995). 

[56] S. Chandrasekhar, Hydrodynamic and Hydromagnetic stability (Clarendon Press, Ox- 
ford, 1961). 

[57] H. Lamb, Hydrodynamics 6th edn (Dover, New York, 1945). 

[58] H. Jeffreys, The Earth, 6th edn (Cambridge University Press, 1976). 

[59] R.R. Lapwood and T. Usami, Free Oscillations of the Earth (Cambridge University 
Press, 1981). 

[60] K. Aki and P.G. Richards, Quantitative Seismology (University Science Books, 
Sausalito, California, 2002). 

[61] J. Hernandez-Roses, M. Picquart, E. Haro-Poniatowski, M. Kanehisa, M. Jouanne 
and J.F. Morhange, Elastic vibrations of spheroidal nanometric particles, J. Phys.: 
Cond. Matter 15 (2003) 7481-7487. 



37 



[62] E. Duval, L. Saviot, A. Mermet and D.B. Murray, Continuum elastic sphere vibra- 
tions as a model for low lying optical modes in icosahedral quasicrystals, J. Phys.: 
Cond. Matter 17 (2005) 3559-3565. 

[63] S.I. Bastrukov and Pik-Yin Lai, Resonant response of a metallic nanoparticle by 
collective cyclotron oscillations of electrons and ions, Phys. Lett. A341 (2005) 207- 
211. 

[64] S.I. Bastrukov, Low frequency elastic response of a spherical particle, Phys. Rev. E49 
(1994) 3166-3170. 

[65] Zhen Ye, On the low frequency elastic response of a spherical particle, Chinese J. 
Phys. 38 (2000) 103-110. 

[66] R. Hasse, Approaches to nuclear friction, Rep. Prog. Phys. 41 (1978) 1027-1101. 

[67] A.L. Fetter and J.D. Walecka, Theoretical Mechanics of Particles and Continua 
(McGraw-Hill, New York, 1980). 

[68] K. Alder and R. M. Steffen, Emission and Absorption of electromagnetic radiation, 
in The Electromagnetic Interaction in Nuclear Spectroscopy (North-Holland, Ams- 
terdam, 1975) 1-55. 

[69] E. Flowers and N. Itoh, Transport properties of dense matter, Astrophys. J. 206 
(1976) 218-242. 

[70] C. Cutler and L. Lindblom, The effect of viscosity on neutron star oscillations, As- 
trophys. J. 314 (1987) 234-241. 

[71] S.I. Bastrukov, J. Yang, D.V. Podgainy and F. Weber, Signatures of field induced 
spin polarization of neutron star matter in seismic vibrations of paramagnetic neutron 
star, J. Phys. G29 (2003) 683-695. 

[72] S.I. Bastrukov, On elastic vibrations of homogeneous star, Mod. Phys. Lett. A8 
(1993) 711-714. 

[73] S.I. Bastrukov, Non-radial vibrations of a star modeled by a heavy spherical mass of 
an elastic solid, Phys. Rev. (1996) E53 1917-1922. 

[74] W. de Boer, C. Sander, V. Zhukov, A. Gladyshev and D. Kazakov, EGRET excess 
of diffuse galactic gamma rays as tracer of dark matter, Astron. and Astrophys. 444 
(2005) 51-67. 



38 



